Forecasting adverse surgical events using self-supervised transfer learning for physiological signals.

Hundreds of millions of surgical procedures take place annually across the world, which generate a prevalent type of electronic health record (EHR) data comprising time series physiological signals. Here, we present a transferable embedding method (i.e., a method to transform time series signals into input features for predictive machine learning models) named PHASE (PHysiologicAl Signal Embeddings) that enables us to more accurately forecast adverse surgical outcomes based on physiological signals. We evaluate PHASE on minute-by-minute EHR data of more than 50,000 surgeries from two operating room (OR) datasets and patient stays in an intensive care unit (ICU) dataset. PHASE outperforms other state-of-the-art approaches, such as long-short term memory networks trained on raw data and gradient boosted trees trained on handcrafted features, in predicting six distinct outcomes: hypoxemia, hypocapnia, hypotension, hypertension, phenylephrine, and epinephrine. In a transfer learning setting where we train embedding models in one dataset then embed signals and predict adverse events in unseen data, PHASE achieves significantly higher prediction accuracy at lower computational cost compared to conventional approaches. Finally, given the importance of understanding models in clinical applications we demonstrate that PHASE is explainable and validate our predictive models using local feature attribution methods.

Autosomal dominant mannose-binding lectin deficiency is associated with worse neurodevelopmental outcomes after cardiac surgery in infants.

OBJECTIVES: The MBL2 gene is the major genetic determinant of mannose-binding lectin (MBL)-an acute phase reactant. Low MBL levels have been associated with adverse outcomes in preterm infants. The MBL2Gly54Asp missense variant causes autosomal dominant MBL deficiency. We tested the hypothesis that MBL2Gly54Asp is associated with worse neurodevelopmental outcomes after cardiac surgery in neonates. METHODS: This is an analysis of a previously described cohort of patients with nonsyndromic congenital heart disease who underwent cardiac surgery with cardiopulmonary bypass before age 6 months (n = 295). Four-year neurodevelopment was assessed in 3 domains: Full-Scale Intellectual Quotient, the Visual Motor Integration development test, and the Child Behavior Checklist to assess behavior problems. The Child Behavior Checklist measured total behavior problems, pervasive developmental problems, and internalizing/externalizing problems. A multivariable linear regression model, adjusting for confounders, was fit. RESULTS: MBL2Gly54Asp was associated with a significantly increased covariate-adjusted pervasive developmental problem score (ß = 3.98; P = .0025). Sensitivity analyses of the interaction between age at first surgery and MBL genotype suggested effect modification for the patients with MBL2Gly54Asp (Pinteraction = .039), with the poorest neurodevelopment outcomes occurring in children who had surgery earlier in life. CONCLUSIONS: We report the novel finding that carriers of MBL2Gly54Asp causing autosomal dominant MBL deficiency have increased childhood pervasive developmental problems after cardiac surgery, independent of other covariates. Sensitivity analyses suggest that this effect may be larger in children who underwent surgery at earlier ages. These data support the role of nonsyndromic genetic variation in determining postsurgical neurodevelopment-related outcomes in children with congenital heart disease.

A vascular endothelial growth factor A genetic variant is associated with improved ventricular function and transplant-free survival after surgery for non-syndromic CHD.

BACKGROUND: We have previously shown that the minor alleles of vascular endothelial growth factor A (VEGFA) single-nucleotide polymorphism rs833069 and superoxide dismutase 2 (SOD2) single-nucleotide polymorphism rs2758331 are both associated with improved transplant-free survival after surgery for CHD in infants, but the underlying mechanisms are unknown. We hypothesised that one or both of these minor alleles are associated with better systemic ventricular function, resulting in improved survival. METHODS: This study is a follow-up analysis of 422 non-syndromic CHD patients who underwent neonatal cardiac surgery with cardiopulmonary bypass. Echocardiographic reports were reviewed. Systemic ventricular function was subjectively categorised as normal, or as mildly, moderately, or severely depressed. The change in function was calculated as the change from the preoperative study to the last available study. Stepwise linear regression, adjusting for covariates, was performed for the outcome of change in ventricular function. Model comparison was performed using Akaike's information criterion. Only variables that improved the model prediction of change in systemic ventricular function were retained in the final model. RESULTS: Genetic and echocardiographic data were available for 335/422 subjects (79%). Of them, 33 (9.9%) developed worse systemic ventricular function during a mean follow-up period of 13.5 years. After covariate adjustment, the presence of the VEGFA minor allele was associated with preserved ventricular function (p=0.011). CONCLUSIONS: These data support the hypothesis that the mechanism by which the VEGFA single-nucleotide polymorphism rs833069 minor allele improves survival may be the preservation of ventricular function. Further studies are needed to validate this genotype-phenotype association and to determine whether this mechanism is related to increased vascular endothelial growth factor production.

Association of Arrhythmia-Related Genetic Variants With Phenotypes Documented in Electronic Medical Records.

IMPORTANCE: Large-scale DNA sequencing identifies incidental rare variants in established Mendelian disease genes, but the frequency of related clinical phenotypes in unselected patient populations is not well established. Phenotype data from electronic medical records (EMRs) may provide a resource to assess the clinical relevance of rare variants. OBJECTIVE: To determine the clinical phenotypes from EMRs for individuals with variants designated as pathogenic by expert review in arrhythmia susceptibility genes. DESIGN, SETTING, AND PARTICIPANTS: This prospective cohort study included 2022 individuals recruited for nonantiarrhythmic drug exposure phenotypes from October 5, 2012, to September 30, 2013, for the Electronic Medical Records and Genomics Network Pharmacogenomics project from 7 US academic medical centers. Variants in SCN5A and KCNH2, disease genes for long QT and Brugada syndromes, were assessed for potential pathogenicity by 3 laboratories with ion channel expertise and by comparison with the ClinVar database. Relevant phenotypes were determined from EMRs, with data available from 2002 (or earlier for some sites) through September 10, 2014. EXPOSURES: One or more variants designated as pathogenic in SCN5A or KCNH2. MAIN OUTCOMES AND MEASURES: Arrhythmia or electrocardiographic (ECG) phenotypes defined by International Classification of Diseases, Ninth Revision (ICD-9) codes, ECG data, and manual EMR review. RESULTS: Among 2022 study participants (median age, 61 years [interquartile range, 56-65 years]; 1118 [55%] female; 1491 [74%] white), a total of 122 rare (minor allele frequency <0.5%) nonsynonymous and splice-site variants in 2 arrhythmia susceptibility genes were identified in 223 individuals (11% of the study cohort). Forty-two variants in 63 participants were designated potentially pathogenic by at least 1 laboratory or ClinVar, with low concordance across laboratories (Cohen κ = 0.26). An ICD-9 code for arrhythmia was found in 11 of 63 (17%) variant carriers vs 264 of 1959 (13%) of those without variants (difference, +4%; 95% CI, -5% to +13%; P = .35). In the 1270 (63%) with ECGs, corrected QT intervals were not different in variant carriers vs those without (median, 429 vs 439 milliseconds; difference, -10 milliseconds; 95% CI, -16 to +3 milliseconds; P = .17). After manual review, 22 of 63 participants (35%) with designated variants had any ECG or arrhythmia phenotype, and only 2 had corrected QT interval longer than 500 milliseconds. CONCLUSIONS AND RELEVANCE: Among laboratories experienced in genetic testing for cardiac arrhythmia disorders, there was low concordance in designating SCN5A and KCNH2 variants as pathogenic. In an unselected population, the putatively pathogenic genetic variants were not associated with an abnormal phenotype. These findings raise questions about the implications of notifying patients of incidental genetic findings.

Burden of potentially pathologic copy number variants is higher in children with isolated congenital heart disease and significantly impairs covariate-adjusted transplant-free survival.

OBJECTIVES: Copy number variants (CNVs) are duplications or deletions of genomic regions. Large CNVs are potentially pathogenic and are overrepresented in children with congenital heart disease (CHD). We sought to determine the frequency of large CNVs in children with isolated CHD, and to evaluate the relationship of these potentially pathogenic CNVs with transplant-free survival. METHODS: These cases are derived from a prospective cohort of patients with nonsyndromic CHD (n = 422) identified before first surgery. Healthy pediatric controls (n = 500) were obtained from the electronic Medical Records and Genetic Epidemiology Network, and CNV frequency was contrasted for CHD cases and controls. CNVs were determined algorithmically; subsequently screened for >95% overlap between 2 methods, size (>300 kb), quality score, overlap with a gene, and novelty (absent from databases of known, benign CNVs); and separately validated by quantitative polymerase chain reaction. Survival likelihoods for cases were calculated using Cox proportional hazards modeling to evaluate the joint effect of CNV burden and known confounders on transplant-free survival. RESULTS: Children with nonsyndromic CHD had a higher burden of potentially pathogenic CNVs compared with pediatric controls (12.1% vs 5.0%; P = .00016). Presence of a CNV was associated with significantly decreased transplant-free survival after surgery (hazard ratio, 3.42; 95% confidence interval, 1.66-7.09; P = .00090) with confounder adjustment. CONCLUSIONS: We confirm that children with isolated CHD have a greater burden of rare/large CNVs. We report a novel finding that these CNVs are associated with an adjusted 2.55-fold increased risk of death or transplant. These data suggest that CNV burden is an important modifier of survival after surgery for CHD.

Next-generation Sequencing of RYR1 and CACNA1S in Malignant Hyperthermia and Exertional Heat Illness.

BACKGROUND: Variants in RYR1 are associated with the majority of cases of malignant hyperthermia (MH), a form of heat illness pharmacogenetically triggered by general anesthetics, and they have also been associated with exertional heat illness (EHI). CACNA1S has also been implicated in MH. The authors applied a targeted next-generation sequencing approach to identify variants in RYR1 and CACNA1S in a cohort of unrelated patients diagnosed with MH susceptibility. They also provide the first comprehensive report of sequencing of these two genes in a cohort of survivors of EHI. METHODS: DNA extracted from blood was genotyped using a "long" polymerase chain reaction technique, with sequencing on the Illumina GAII or MiSeq platforms (Illumina Inc., USA). Variants were assessed for pathogenicity using bioinformatic approaches. For further follow-up, DNA from additional family members and up to 211 MH normal and 556 MH-susceptible unrelated individuals was tested. RESULTS: In 29 MH patients, the authors identified three pathogenic and four novel RYR1 variants, with a further five RYR1 variants previously reported in association with MH. Three novel RYR1 variants were found in the EHI cohort (n = 28) along with two more previously reported in association with MH. Two other variants were reported previously associated with centronuclear myopathy. The authors found one and three rare variants of unknown significance in CACNA1S in the MH and EHI cohorts, respectively. CONCLUSIONS: Targeted next-generation sequencing proved efficient at identifying diagnostically useful and potentially implicated variants in RYR1 and CACNA1S in MH and EHI.

Patient genotypes impact survival after surgery for isolated congenital heart disease.

BACKGROUND: Survival after cardiac surgery in infancy requires adaptive responses from oxidative stress management and vascular regulation pathways. We tested the hypothesis that genetic variation in these pathways influences postoperative survival in nonsyndromic congenital heart disease children. METHODS: This is an analysis of a cohort of nonsyndromic congenital heart disease patients who underwent cardiac surgery with cardiopulmonary bypass before 6 months of age (n=422). Six single nucleotide polymorphisms (SNPs) in six genes involved in oxidative stress and vascular response pathways, identified through a priori literature search, were tested for effects on transplant-free survival. Survival curves, adjusting for confounding covariates, were calculated using the Cox proportional hazard models. RESULTS: Long-term survival was strongly associated with vascular endothelial growth factor A gene SNP rs833069 (p=7.03×10(-4)) and superoxide dismutase 2 gene SNP rs2758331 (p=0.019). To test for joint effects of the two SNPs on transplant-free survival, the genotypes were grouped to form a risk score reflecting the cumulative number of risk alleles (0 to 4 alleles per patient). A higher risk score based on the VEGFA and SOD2 SNP genotypes was associated with worse transplant-free survival (p=3.02×10(-4)) after confounder adjustment. The total burden of risk alleles was additive; subjects with the highest risk score of 4 (n=59 subjects, 14.2% of the cohort) had a total covariate-adjusted hazard ratio of 15.64 for worse transplant-free survival. CONCLUSIONS: After cardiac surgery, infants who are homozygous for the high-risk alleles for both the VEGFA and SOD2 SNPs have an approximately 16-fold increased risk of death or heart transplant, suggesting that genetic variants are important modifiers of survival after surgery for congenital heart disease.

Actionable, pathogenic incidental findings in 1,000 participants' exomes.

The incorporation of genomics into medicine is stimulating interest on the return of incidental findings (IFs) from exome and genome sequencing. However, no large-scale study has yet estimated the number of expected actionable findings per individual; therefore, we classified actionable pathogenic single-nucleotide variants in 500 European- and 500 African-descent participants randomly selected from the National Heart, Lung, and Blood Institute Exome Sequencing Project. The 1,000 individuals were screened for variants in 114 genes selected by an expert panel for their association with medically actionable genetic conditions possibly undiagnosed in adults. Among the 1,000 participants, 585 instances of 239 unique variants were identified as disease causing in the Human Gene Mutation Database (HGMD). The primary literature supporting the variants' pathogenicity was reviewed. Of the identified IFs, only 16 unique autosomal-dominant variants in 17 individuals were assessed to be pathogenic or likely pathogenic, and one participant had two pathogenic variants for an autosomal-recessive disease. Furthermore, one pathogenic and four likely pathogenic variants not listed as disease causing in HGMD were identified. These data can provide an estimate of the frequency (∼3.4% for European descent and ∼1.2% for African descent) of the high-penetrance actionable pathogenic or likely pathogenic variants in adults. The 23 participants with pathogenic or likely pathogenic variants were disproportionately of European (17) versus African (6) descent. The process of classifying these variants underscores the need for a more comprehensive and diverse centralized resource to provide curated information on pathogenicity for clinical use to minimize health disparities in genomic medicine.

Exome sequencing reveals novel rare variants in the ryanodine receptor and calcium channel genes in malignant hyperthermia families.

BACKGROUND: About half of malignant hyperthermia (MH) cases are associated with skeletal muscle ryanodine receptor 1 (RYR1) and calcium channel, voltage-dependent, L type, α1S subunit (CACNA1S) gene mutations, leaving many with an unknown cause. The authors chose to apply a sequencing approach to uncover causal variants in unknown cases. Sequencing the exome, the protein-coding region of the genome, has power at low sample sizes and identified the cause of over a dozen Mendelian disorders. METHODS: The authors considered four families with multiple MH cases lacking mutations in RYR1 and CACNA1S by Sanger sequencing of complementary DNA. Exome sequencing in two affecteds per family, chosen for maximum genetic distance, were compared. Variants were ranked by allele frequency, protein change, and measures of conservation among mammals to assess likelihood of causation. Finally, putative pathogenic mutations were genotyped in other family members to verify cosegregation with MH. RESULTS: Exome sequencing revealed one rare RYR1 nonsynonymous variant in each of three families (Asp1056His, Val2627Met, Val4234Leu), and one CACNA1S variant (Thr1009Lys) in the fourth family. These were not seen in variant databases or in our control population sample of 5,379 exomes. Follow-up sequencing in other family members verified cosegregation of alleles with MH. CONCLUSIONS: The authors found that using both exome sequencing and allele frequency data from large sequencing efforts may aid genetic diagnosis of MH. In a sample selected by the authors, this technique was more sensitive for variant detection in known genes than Sanger sequencing of complementary DNA, and allows for the possibility of novel gene discovery.

Identifying safer anti-wear triaryl phosphate additives for jet engine lubricants.

Individuals aboard jet aircraft may be exposed to potentially toxic triaryl organophosphate anti-wear lubricant additives (TAPs) that are converted by cytochromes P450 into toxic metabolites. Consequences of exposure could be reduced by using less toxic TAPs. Our goal was to determine whether an in vitro assay for inhibition of butyrylcholinesterase (BChE) by bioactivated TAPs would be predictive of inhibition of serine active-site enzymes in vivo. The in vitro assay involved TAP bioactivation with liver microsomes and NADPH, followed by incubation with human BChE and measurement of BChE activity. Of 19 TAPs tested, tert-butylated isomers produced the least BChE inhibition. To determine the relevance of these results in vivo, mice were exposed to Durad 125 (D125; a commercial mixture of TAP esters) or to TAPs demonstrating low or no BChE inhibition when assayed in vitro. Inhibition of BChE by bioactivated TAPs in vitro correlated well with inhibition of other serine active-site enzymes in vivo, with the exception of brain acetylcholinesterase and neuropathy target esterase (NTE), which were not inhibited by any TAP tested following single exposures. A recombinant catalytic domain of NTE (rNEST) exhibited classical kinetic properties of NTE. The metabolite of tri-(o-cresyl) phosphate (ToCP), 2-(o-cresyl)-4H-1,3,2-benzodioxaphosphoran-2-one (CBDP), inhibited rNEST in vitro, but with an IC(50) value almost 6-times higher than for inhibition of BChE. Physiologically-relevant concentrations of the flavonoid naringenin dramatically reduced D125 bioconversion in vitro. The in vitro assay should provide a valuable tool for prescreening candidate TAP anti-wear additives, identifying safer additives and reducing the number of animals required for in vivo toxicity testing.

Biomarkers of organophosphorus (OP) exposures in humans.

There are ongoing events where aircraft engine lubricant containing tricresyl phosphates (TCPs) contaminates aircraft cabins. Some individuals have experienced tremors or other neurological symptoms that may last for many months following exposures. Mass spectrometric (MS) protocols are being developed to determine the percentage of "biomarker proteins" that are modified by such exposures, specifically on active site serines. Both plasma butyrylcholinesterase (BChE) and red cell acylpeptide hydrolase (APH) are readily inhibited by 2-(ortho-cresyl)-4H-1,3,2-benzodioxaphosphoran-2-one (CBDP) or phenyl saligenin cyclic phosphate (PSP) and have the potential to provide information about the level of exposure of an individual. We have developed immunomagnetic bead-based single-step purification protocols for both BChE and APH and have characterized the active site serine adducts of BChE by MS.

Identification and characterization of biomarkers of organophosphorus exposures in humans.

Over 1 billion pounds of organophosphorus (OP) chemicals are manufactured worldwide each year, including 70 million pounds of pesticides sprayed in the US. Current methods to monitor environmental and occupational exposures to OPs such as chlorpyrifos (CPS) have limitations, including low specificity and sensitivity, and short time windows for detection. Biomarkers for the OP tricresyl phosphate (TCP), which can contaminate bleed air from jet engines and cause an occupational exposure of commercial airline pilots, crewmembers and passengers, have not been identified. The aim of our work has been to identify, purify, and characterize new biomarkers of OP exposure. Butyrylcholinesterase (BChE) inhibition has been a standard for monitoring OP exposure. By identifying and characterizing molecular biomarkers with longer half-lives, we should be able to clinically detect TCP and OP insecticide exposure after longer durations of time than are currently possible. Acylpeptide hydrolase (APH) is a red blood cell (RBC) cytosolic serine proteinase that removes N-acetylated amino acids from peptides and cleaves oxidized proteins. Due to its properties, it is an excellent candidate for a biomarker of exposure. We have been able to purify APH and detect inhibition by both CPS and metabolites of TCP. The 120-day lifetime of the RBC offers a much longer window for detecting exposure. The OP-modified serine conjugate in the active site tryptic peptide has been characterized by mass spectrometry. This research uses functional proteomics and enzyme activities to identify and characterize useful biomarkers of neurotoxic environmental and occupational OP exposures.